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Since the discovery that the copper oxide sheets in the High Tc superconductors show- 
strong antiferromagnetic (AF) spin correlations [?], there has been increased interest in 
studying magnetic properties of these systems. This followed from the argument, made by 
Anderson and many other researchers [?], that the essential physics of the copper oxide 
superconductors is contained in doped antiferromagnets, where a central issue is the rela- 
tionship between the hole doping and AF spin correlations [?,?]. In particular, the phase 
diagrams of these materials as functions of hole doping have been established by neutron 
scattering, muon spin rotation, and magnetic resonance measurements [?]. When the hole 
doping concentration exceeds some critical value (about 5%), AF long-range order (AFLRO) 
disappears and the materials are converted into nonmagnetic metals. It is believed [?] that 
the physics of these materials may be effectively described by a 2D, large U Hubbard model 
or its equivalent, the t-J model. 

In an attempt to understand the relationship between the hole doping and AF spin cor- 
relations, many authors have studied holes moving in the background of the spin resonating 
valence bond state [?], spin flux phase [?], and spiral spin phase [?]. A hole in the AF 
background has a very large effective mass because of frustrations [?]. Numerical simula- 
tions suggest [?] that the 2D t-J model is Neel ordered at half-fiUing at temperature T=0, 
and the suppression of AFLRO as a consequence of hole doping has already been studied 
by using variational Monte-Carlo technique [?] . Recently, the fermion-spin theory has been 
employed to study the ground-state properties of the 2D t-J model [?,?] within the mean- 
field-approximation (MFA). The results show that the magnetized vr-flux state, which is the 
coexistence of AFLRO with a vr-flux state, has the lowest energy at half-fiUing, and the 
AFLRO is destroyed by hole doping of the order ~ 10% ~ 15% for reasonable values of 
the parameters. However, quantum fluctuations are dropped completely in this mean-field 
treatment. In fact, there is a strong coupling between spin (spinon) and charge (holon) de- 
grees of freedom [?] . It has been emphasized [?] that for small spin values, strong quantum 
fluctuations of spinous may generate a novel spin liquid ground-state. The purpose of this 
paper is to study these issues within the t-J model by considering quantum fluctuations of 
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spinons. 

We start from the t-J model on a square lattice which can be written as 

H — —t ^ C/^Cj+^CT + h.c. — // ^ Cl^Cia + J ^ Sj • Si+^, (1) 

icr,T]>0 ia i,ri>0 

with r] — x,y, C}^ {Cia) are the electron creation (annihilation) operators, Sj = CjaCi/2 
are spin operators with a — {ax, Cy, Uz) as Pauli matrices, and // is the chemical potential. 
The t-J Hamiltonian (1) is supplemented by the on-site local constraint, Yja^laCia < 1, 
i.e., there should be no doubly occupied sites. 

The crucial requirement for the t-J model is to impose this local constraint for a proper 
understanding of the physical properties of the copper oxide materials [?] . The local nature 
of the constraint is of prime importance, and its violation may lead to unphysical results [?] . 
To avoid these difficulties, a fermion-spin transformation based on the charge-spin separation 

= h\S., Q = h\St, (2) 

was proposed [?] to study the t-J model, where the spinless fermion operator hi keeps track 
of the charge (holon) while the pseudospin operator Si keeps track of the spin (spinon) . In 
this case, the on-site local constraint for single occupancy is satisfied even in the MFA. For 
discussing quantum spin problems within the framework of fermion many particle theory, we 
map the quantum spin onto the spinless fermion representation via the 2D Jordan- Wigner 
[?] transformation, Sf — //e*^% S^' — fie~^^\ with /j as a spinless fermion, phase factor 
— Y^i^i fifiBii, and Bn — Imln(Z/ — Zj), in the complex representation of the lattice sites 
Zi — Xi + iYi. In this case, it has been shown that the t-J model can be expressed as [?,?,?] 

H^Ho + H?^ + hP, (3a) 

i ^ i,T]>0 

-2JeffY.flf^ + \NJef^, (3b) 

H? ^-tY: h^hl^{e'C'''■^'f}h^, + e-^/r A-'^Vi./O + h.c (3c) 
= ^ E (2Je//)///./i,/.+„ (3d) 
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where N is the number of sites, and Jeff ~ J{(hihl)(hi+rjhl_^_^)), H^P and Hf^ are the 
spinon-holon and spinon-spinon interaction Hamiltonians, respectively. A(rj) is the vector 
potential created by the 2D Jordan- Wigner phase factor, which relates the 2D quantum spin 
problem to the fermion's flux phase problem [?]. The Hamiltonian (3) contains both phase 
and amplitude fluctuations. The phase fluctuations can be studied by using the anyon and 
eikonal expansion techniques [?]. However, the essential ingredients to give a qualitatively 
correct description of the global ground-state properties perhaps are included in amplitude 
fluctuations [?]. Therefore, in this paper, we only discuss the amplitude fluctuations and 
treat the phase factor in eq. (3) in MFA, as in the previous works [?,?,?]. For convenience, 
we choose the gauge field 

= (4) 

which is produced by a | flux quantum (or phase tt) through each plaquette. 

For discussing the self-consistent spinon and holon Green's functions, we define the 
nearest-neighbor spin bond-order amplitude, Neel order parameter, and holon particle-hole 
parameter as 

(5f ) = (-l)^M, (6) 
0=(/iJ/i,+,), (7) 

respectively. This means we consider the coexistence of AFLRO with a 7r-flux state, and this 
state, referred to as the magnetized 7r-fiux state, does not break the time-reversal symmetry 
and parity [?,?,?]. The site subscripts of the order parameters x, 0, and M have been 
dropped since the system is translation invariant. 

Since there are two inequivalent sublattices A and B in the AF system, the one-particle 
spinon and holon Matsubara Green's functions can be defined as matrices, 

D(k, t-t') = Di(k, T - r') + Dt^K t-t'), (8a) 
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{Trfi(r)fj^\r')) 

(7;/f(r)/f(r')) 



(8b) 
(8c) 



and 



^(k, T - r') = ^^(k, T - r') + G'T(k, r - r'), (9a) 
/(r./i^(r)/if(r')) \ 

( (T./.^(r)/.f(r'))\ ^ ^ 

Gt k, T - r' = - , 9c 

where r and r' are the imaginary times, and (. . .) is an average over the ensemble. Dl{Gl) 
and Dt{Gt) are the longitudinal and transverse parts of the spinon (holon) Green's function 
D{G), respectively. 

A self-consistent method of treating interaction Hamiltonians as perturbation has been 
given in Ref. [?] within the Schwinger-boson slave-fermion theory. We [?] have used similar 
self-consistent calculations within the Jordan- Wigner approach to discuss the t-J model at 
half-filling and obtained some interesting results. In the following, we will generalize the 
method in Ref. [?] to discuss the t-J model away from half- filling within the fermion-spin 
approach. 

With the gauge choice of eq. (4), the zero-order spinon and holon Green's functions of 
Hq in eq. (3b) are given by 

- {iu^ + 2 J.,,)^ - u;g(k) ' - {iu^ + 2 Je,,)^ - a;o^(k) ' ^'"^ 

and 

Gf{k)^-^, G^^\k)^0, (11) 

respectively, where k — (k, ia;„), a;o(k) = Jey/Fk, and Fk = ^Jcos^k^ -\- cos^/c^. 

According to the Dyson's theory [?], the full spinon and holon Green's functions of H 
in eq. (3) due to considering the spinon-holon and spinon-spinon interactions are formally 
expressed as. 



^^^^ ^(o)-i(-/j) _ j:^[k) iujn + 2Je// - Jefficr^cosk^ - aycosky) - T^s{k) ' ^""^^^ 

^^^^ " GW-HA;) - E;,(A;) " + - E;,(A;) ' ^^^^ 

where ^^(A;) and S/i(A;) are the spinon and holon self-energy functions, respectively. In the 
self-consistent calculation, the spinon and holon Green's functions cannot be self-consistently 
determined only by Eqs. (12) and (13), since the order parameters x, M, and (p are also 
included in the spinon and holon self-energy functions, and we must have another set of 
self-consistent equations to determine these parameters. According to the definitions of the 
order parameters in Eqs. (5)- (7) and the gauge choice of (4), we can get the following 
equations, 
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X = -gY.Tr[{a^cosk^ - ayCOsky)DT{'k, iuj^)], (14) 

M=~ i: WTr[a,DL(\^,iu;n)], (15) 
<^ = ^1 E l^^j:Tr[a,GTik,iu;^)], (16) 
^=11 E ^T.TrG,{k,iu;„), (17) 

^ -'^ k(r-ed) ^ iu>n 

where (red) means the summation is carried only over the reduced Brillouin Zone. The order 

parameters M, (p, and chemical potential jj, can be obtained by self-consistently solving 
Eqs. (14)-(17). 

Although the matrix analysis is quite tedious, the calculation of the first-order self-energy 

(Hartrcc-Fock self-energy) function is straightforward. Following the discussions in Ref. [?], 
we obtain the spinon Hartree-Fock self-energy 

^iL^\k) = 2Je// - mpa^, E''J^^\k) = 2{(j)t - xJeff){crxCosk^ - ayCOsky), (18) 

and the holon Hartree-Fock self-energy 

^iT\k) = 0, E(?^)(A;) = ek<7„ (19) 



where p = (1 - 2x)Je// + 2t0, m = IM^^ ^ 4xi7k, and 7k = KcosA;^; + cosky). 

) and S^^^^(e|^^^'*) are the longitudinal and transverse parts of the spinon 
(holon) Hartree-Fock self-energy E^^^)(S^^^''), respectively. Substituting these results of 
^iL^\k), ^iT^\k), and T.f/\k) into Eqs. (12) and (13), we obtain self- 

consistent Hartree-Fock (SCHF) spinon and holon Green's functions, 

D^"^\k) = ^^»-mp^. D^^^\k) = - ^,,cosfc^) ^ 

and 



respectively, where a;(k) = pE'k and £'k = -y/F^ -|- m^. This SCHF solution is exactly the 
usual self-consistent MFA in the previous works [?,?]. At half-filling, the staggered magne- 
tization and ground-state energy are M — 0.389 and Eq — — 0.33J per bond, respevtively. 
Away from half-fiUing, AFLRO vanishes around doping 5 = 10 — 15% for reasonable values 
of the parameter t/J, which are in rough agreement with experiments [?] and Monte- Carlo 
simulations [?] . In this case, the SCHF Hamiltonian describes free spinous and holons mov- 
ing in an average Chern-Simons gauge field of flux tt per plaquette [?]. 

The mean-field phase boundary between AFLRO and disordered states is, of course, at 
somewhat higher doping S than the value given by experiments and numerical simulations. 
The t-J model is characterized by a competition between the kinetic energy (t) and magnetic 
energy (J). The magnetic energy J favors an AFLRO for spins, whereas the kinetic energy 
t favors delocalization of holes and tends to destroy AFLRO. According to the experiments 
[?], the suppression of AFLRO depends on both hole dopings and quantum fluctuations of 
spinous. Therefore in this paper the holon part will be restricted to the Hartree-Fock level, 
and we consider quantum fluctuations of spinous due to the spinon-spinon interaction Hj^^ 
in Eq. (3d) within the random-phase- approximation (RPA). 

The analytic spinon Green's function is very complicated as emphasized in Refs. [?,?]. 
To make the actual calculation feasible, we start from the SCHF solution, i.e., the full spinon 



Green's function is replaced by the spinon Hartree-Fock Green's function in the calculation 
of the spinon RPA self-energy [?,?]. The detailed process to obtain the spinon RPA Green's 
function within the t-J model at half-filling has been given in Ref. [?]. Following the pocedure 
described there [?], we get the spinon RPA Green's function in the present case, 
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iujn — [mp — A(k)]az + [p + B(k)]{aj;COskx — Oycosky) 



(22) 



where ujRPA{k) = ^[mp - A{k)]^ + [p + B{k)]^ri, and 

-{uj^ + a;p+q + a;k+p) + ^^k 



X — — ^ — , (23a) 

(a;q u;p+q + a;k+p)^ + oj^ 



(a;q + a;p^q + a;k+ p)+c.k 



Here Srpa{p) is the spinon RPA dielectric function. With the above holon and spinon 
Green's functions, we obtain the self-consistent equations at zero temperature , 

rnp-Aik^ (25) 
</'=^ E ^[e{p-e^)-e{p + e^)], (26) 

k(rec() 

^ = 1^ E ^[^(/^-ek) + ^(/^ + 6k)], (27) 

k(re(i) 

where the effect of quantum fluctuations of spinous from the spinon-spinon interaction have 
been included in A{k), B{k), and u!RPA{k), which will affect the parameters x, M, and p, 
but not the parameter at low dopings. This is because the effect due to the variance of 
parameter x in ek = 4xt7k in Eqs. (26) and (27) can be incorporated into the parameter t/J. 
On the other hand, we have shown in Ref. [?] that at low dopings the holon particle-hole 
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order parameter increases roughly linearly upon doping, and is almost independent oit/ J. 
Therefore in this case we can take the mean-field values of the parameter cj) to calculate the 
RPA self-consistent equations (24) and (25) for the parameters x and M. We have performed 
such RPA calculation. At half-fiUing, the result is the same as those in Ref. [?], i.e. the 
staggered magnetization and RPA ground-state energy are M — 0.327 and Eq — — 0.332J 
per bond, respectively, which are very close to the best numerical result [?] M = 0.31 and 
Eo — — 0.3346J. Away from half-filling, AFLRO is quickly destroyed by doping of the order 
oi 5 — 5% for a reasonable parameter value t/J — 5, which is plotted in Fig. 1 (solid 
triangles). For comparison, the experimental measurements on {Lai-xBax)2Cu04 (solid 
squares, the dashed fine is guide for eyes) [?] and numerical results (sohd circles) [?] are 
also plotted in Fig. 1. From Fig. 1, it is obvious that our result is in agreement with the 
experiments [?] and Monte-Carlo simulations [?]. As compared with the mean- field result 
[?], the present RPA results indicate that both hole dopings and quantum fiuctuations of 
spinous lead to a strong suppression of AFLRO. 

In summary, we have studied the doping dependence of AFLRO in the copper oxdie 
materials within the t-J model. We have performed a self-consistent calculation, where the 
holon part is restricted to the Hartree-Fock level, and the quantum fiuctuation of spinous is 
considered within the RPA. We have considered only the amplitude fiuctuations. AFLRO 
vanishes around doping S — 5% for a reasonable parameter value t/J — 5, which is in 
agreement with the experiments on copper oxide materials [?] and Monte-Carlo simulations 

[?]• 
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